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I. INTRODUCTION 

One of the major goals of condensed-matter physics 
during the past decades has been to understand the role 
of electronic correlations in solids. At the heart of this 
field lies the metal-insulator transition (MIT) driven by 
the electronic interaction, the Mott transition^ It occurs 
in many materials^ in particular transition-metal com- 
pounds such as V2O3. The underlying physical mech- 
anism for this transition is captured by the single-band 
Hubbard model; 3 - 
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Here cJ CT are creation operators for an electron at site 
i with spin a, U is the local Coulomb repulsion, and 
Uj arc the hopping amplitudes, e.g., t\ for hopping be- 
tween nearest neighbors (NN) and ti for hopping be- 
tween next-nearest neighbors (NNN). In the Hubbard 
model, the Mott transition can be driven either by an 
increase of the interaction U (at half- filling) , or, for 
large enough U, by a change of the particle density to- 
wards one electron per site. The two cases are referred 
to as bandwidth-controlled Mott transition and filling- 
controlled Mott transition, respectively^ 

Much progress in understanding electronic correla- 
tion effects in general, and the Mott transition in par- 
ticular, was made by considering the limit of infinite 
dimensions^ Dynamical mean-field theory (DMFT)££ 
can then be used as a unified framework to study metal- 
lic, insulating, and magnetically ordered phases. Both 
the bandwidth-controlled MIT 5 i 7 i 8 iV Q i n i 12 i 13 i 14 i 15 and 
the filling-controlled MIT 5 i 16 i 17 i 18 i 19 i 20 i 21 were studied 
within DMFT; for T > the transition is first-order. 
However, in the case of a bipartite lattice and if Uj is re- 
stricted to NN hopping, i.e., in the particle- hole symmet- 
ric case, the Mott transition at half-filling quite generally 
does not occur since it is preempted by the intrinsic weak- 
coupling instability towards antiferromagnetic (AF) long- 
range order. To reduce or even avoid this "parasitic" 



AF phase one may consider non-bipartite lattices 2 ^ 3 - or 
alternatively introduce NNN hopping. The latter in- 
troduces a competition which frustrates the antiferro- 
magnetic order. It was proposed that the phase dia- 
gram of the single-band Hubbard model with NN and 
NNN hopping then qualitatively resembles that of vari- 
ous correlated-electron materials and thus appears to be 
the minimal model for such materials 

In addition to the frustration of the antiferromag- 
netic phase, NNN hopping generally makes the Hub- 
bard model particle-hole asymmetric, even at half-filling. 
This asymmetry, which is a generic property of real 
materialsj2£i 2 & will also have an effect on the paramag- 
netic phase. However, previous DMFT investigations 
considered a special type of hopping ("random <2 hop- 
ping" ) in which the model remains particle- hole symmet- 
ric even for t 2 ^ 0, and the paramagnetic phase is not 
changed at a \\£ i 2 i 21i2L2& By contrast, the effect of non- 
random NNN hopping becomes evident already in the 
non-interacting system through an asymmetric density 
of states (DOS), as derived recently by a new theoreti- 
cal technique which is able to treat arbitrary hopping on 
the Bethe latticei 29 i 30 In the present paper we proceed 
to discuss the effects of NNN hopping t 2 on the phase 
diagram for finite U and T. 

As we will show, a striking consequence of NNN hop- 
ping is the occurrence of phase separation between metal- 
lic and insulating phases even at half-filling. Indeed, in 
the paramagnetic phase of the Hubbard model, phase 
separation is known to lead to hysteresis in the den- 
sity n{n) at the Mott transitionii2i2i For the Hubbard 
model on a square lattice (d = 2), which is closely related 
to high-T c superconductors, phase separation was found 
within the dynamical cluster approximation (DCA) for 
<2 > and n < 1^ whereas for t 2 = finite-size Monte 
Carlo calculations 3 ^ 3 , provide evidence for a homoge- 
neous state. On the other hand, in infinite spatial di- 
mensions DMFT predicts phase separation already for 
t 2 = and n ^ but not at half-filling. We find 

that this changes for ti-t 2 hopping. 
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For our investigation of the effect of ti hopping on 
the paramagnetic phase diagram we use the recently 
developed self-energy functional approach (SEA) ^4 The 
SFA is a variational method based on an exact varia- 
tional principle for the grand potential^ In this scheme, 
the self-energy of a finite reference system is taken as 
variational ansatz for the self-energy of the lattice sys- 
tem. The method is very general and can be applied 
to finite-dimensional systems^ and for various types of 
interactions^ Within DMFT the exact result would be 
recovered for a single-impurity Anderson model (SIAM) 
with infinitely many bath sites as reference system^ 
while in practice only a SIAM with a finite number n s of 
bath sites can be investigated. However, the results con- 
verge quickly with increasing n s £& The simplest impurity 
model considers only a single bath site. Surprisingly, this 
reference system is already enough to reproduce previous 
results for = for the Mott transition qualitatively 
and, sometimes, even quantitatively^ Therefore we will 
use this "two-site approximation" (one correlated site, 
one bath site) also to investigate the case t% ^ 0. 

The paper is organized as follows: In Scc.|H]we gi ye 
a short description of the SFA and of our method to 
treat NNN hopping on the Bethe lattice. In Sec.lTTTlNN 
hopping is studied at and away from half-filling. We cal- 
culate the paramagnetic phase diagram as a function of 
interaction, temperature, and either chemical potential 
or density. The phase diagrams for NN and NNN hop- 
ping arc discussed in detail in Sec. IIVI 



II. MODELS AND METHODS 

A. Self-energy functional approach 

The self-energy functional approach (SFA) is based on 
the variational principle of Luttinger and Wardi^ When 
formulated in terms of the self-energy S, the basic quan- 
tity in this context is the grand-potential functional f2[S], 
which becomes stationary at the physical self-energyiM It 
can be written as 

fi[S] = F[S] + Trlog[-(Go 1 - S)" 1 ] , (2) 

where F[S] contains all information about the interac- 
tion, while the noninteracting Green function Go is con- 
sidered fixed. (Quantities written in boldface are matri- 
ces in single-particle indices and Matsubara frequencies, 
and Tr denotes the corresponding trace.) The functional 
form of F[S] is generally unknown, apart from the fact 
that it depends only on the interaction part in the Hamil- 
tonian, and not on the kinetic energy. This universal- 
ity follows directly from various ways in which F[S] or 
its Legendre transform, the Luttinger- Ward functional, 
can be constructed 4& It allows one to calculate 0[S] ex- 
actly for certain variational self-energies, namely for all S 
that can be considered as the exact self-energy of a refer- 
ence system with the same interaction. The self-energy, 



£ = E(yj) is varied by varying single-particle parame- 
ters {yi} of that reference system. In general X can only 
be calculated for small systems, by means of exact diag- 
onalization. Comparing with the self-energy functional 
of the reference system and exploiting the universality of 
F[S], Eq. J2Jl can be written as: 

fi[S] = £1' - Trlog[— (G _1 - S) -1 ] 

+ Trlog[-(G - 1 -E)- 1 ], (3) 

where VL' and G correspond to the exact grand potential 
and the noninteracting Green function of the reference 
system, respectively. The best approximation for the self- 
energy is then determined from dCl[S(yi)]/dyi = 0. 

As noted in the Introduction, the reference system em- 
ployed in our investigation is a two-site impurity model^ 

H re f = Un^mi+ ^ e j"iCT+^^(q CT C2 CT +h.c.) . (4) 

i— 1,2<7 a 

Here index 1 refers to the correlated impurity and index 
2 to the bath site. The on-site energies ei,2 and the hy- 
bridization V are taken as variational parameters. As we 
are interested only in the paramagnetic phase and thus 
need to consider only spin-independent self-energies, we 
can take and V to be spin-independent as well. 

Once a stationary point of £1 is found, it is tracked in 
the space of variational parameters {yi} as a function of 
the external physical parameters (interaction U, temper- 
ature T, and chemical potential fi) using local algorithms 
for the solution of the equations dSl/dyi = 0. High ac- 
curacy is achieved by calculating derivatives of d^l/diji 
analytically^ However, close to a first-order transition 
it is inconvenient to parametrize the solution by the phys- 
ical parameters, because in this parametrization several 
solutions coexist which correspond to various metastable 
phases. To avoid repeated switching between these so- 
lutions we use an algorithm that treats the physical and 
variational parameters on equal footing. 

B. NN and NNN hopping on the Bethe lattice 

The properties of the lattice enter Eq. only via the 
term Tr logHG^ 1 - I!)" 1 ] = f) laU [£]. In this subsec- 
tion we show how to evaluate Oi a tt f° r arbitrary hop- 
ping on the Bethe lattice, in particular for t\-t2 hopping. 
We restrict the analysis to the paramagnetic phase with 
a spin- and site-independent self-energy S = T,{iuj n ). 
Generally the trace in f2i a tt includes a sum over Matsub- 
ara frequencies and over all eigenstates \k) of the matrix 
Gq 1 — £ = iu n + 1-1 — £(wj n ) — t, where t is the matrix of 
hopping amplitudes. The sum over \k) is then expressed 
in terms of a sum over lattice sites i and an integral over 
the local density of states, pj(e) = J2 k S(e — ek)\(i\k)\ 2 , 
and the sum over Matsubara frequencies is transformed 
into an integral over real frequencies as done in Ref. I381 
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This leads to the expression 

oo oo 

^latt = 2 J J de pi(e)Q[uj + fx - - e] , 

1 — OO —OO 

(5) 

where f(u>) is the Fermi function and the factor 2 ac- 
counts for spin degeneracy. For the Bethe lattice with 
arbitrary coordination number Z there exists a general 
method to calculate the DOS for any given hopping 
Hamiltoniani^iSi One first introduces hopping Hamil- 
tonians Hd = Y]j. ,_j that describe hopping only 
between dth nearest neighbors. In the limit of infinite 
connectivity K = Z — 1 these Hamiltonians must be 
scaled^ according to Hd = Hd/K d / 2 , where Hd retains 
a nontrivial spectrum for K — *■ oo. In particular, the 
spectrum of the NN hopping Hamiltonian Hi has the 
well-known semi-elliptical form po{\) = V4 — A 2 / (2-k) in 
this limit. One may then use a remarkable relation valid 
for the Bethe lattice, which makes use of its special topo- 
logical properties^ Namely, every hopping Hamiltonian 
can be written as a function of Hi, i.e., for an arbitrary 
set of hopping amplitudes {t d } one has 

Y,m d = e{H l ). (6) 

d 

Analytical expressions for e(A) were derived in Ref. I29I 
In analogy to crystal lattices, where any translationally 
invariant Hamiltonian is a function of momentum, e(A) 
plays the role of a dispersion relation, and A € [—2, 2] 
runs over the spectrum of Hi. In particular, the disper- 
sion for t\-t2 hopping is given by2S 

e(A)=t;A 2 + ^A-t;. (7) 

The DOS can then be obtained in a straightforward way 
from the known DOS for Hi by a change of variables. 
The lattice contribution for arbitrary range hopping then 
reduces to 



niatt= / dujf(tu) I d\^-— e[w+/i-E(w)-e t .(A)]. 

-oo -2 

(8) 

This expression is suitable for numerical evaluation. For 
t\-t2 hopping the inner integral can be evaluated analyt- 
ically by solving a quadratic equation. This allows us to 
obtain Oi a tt with great precision, which is necessary to 
determine the stationary points of £1 reliably. 

We conclude this section with a general remark about 
the sign of the hopping amplitudes t\ and £2- For any bi- 
partite lattice the unitary transformation Cj CT — > (— l) 1 ^, 
where (— 1)* is alternating on the two sublattices, changes 
the sign of all hopping amplitudes between the two sub- 
lattices. Thus we can assume ti > without loss of gen- 
erality. On the other hand, it follows from the particle- 
hole transformation Cj CT — > (— l) l c| CT that <2 — ► —^2 merely 
interchanges electron-doped (n > 1) and hole-doped 
(n < 1) regimes. Therefore we can also assume ti > 0. 
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FIG. 1: (U,fi) phase diagram at T = for NN hopping. Re- 
sults from two-site SFA (this work) are compared with ED 
(EDI: 8 sitesii, ED2: extrapolation from 11 sitea^) and 
NRG^- Due to particle-hole symmetry the phase diagram 
is symmetric with respect to /j, — U /2. 



III. NN HOPPING 

The paramagnetic phase-diagram of the Hubbard 
model with only NN hopping (£2 = 0) has been stud- 
ied intensively using DMFT, both for half-nllingVV 3 ! 14 
and for finite dopingAiSiiLiSiSSiSi In this section we use 
the two-site approximation to investigate the full param- 
eter space spanned by U, T, and /i, and compare with 
previous results where available. 

At T = the phase diagram has a simple structure 
(Fig. UJ). At half-filling (where fx = U/2) the system 
is insulating for U > U C 2, where U C 2 ~ 5.84t^ for the 
semi-elliptic DOS>i£ For any other filling the ground 
state is metallic. In the metal the charge compressibility 
Xfi = dn/dfi remains finite even when the insulator is 
approached from larger and lower filling, but the chem- 
ical potential changes discontinuously from /i = fj,~ 2 to 
/i = as the filling goes from n = 1~ to n = 1 + . For 
fx~ 2 < M < Mc2 the system is insulating and half-filled 
(Fig. P). In a region < // < /i~ 2 an( l M^i > M > Mc2 
both metallic and insulating DMFT solutions exist, with 
the metallic phase being the thermodynamically stable 
one. The appearance of the metal at /1 = /i^ lii related 
to the development of states inside the Mott gap, whereas 
the breakdown of the insulator at /1 = fi cl occurs when 
/i reaches the edge of the gap^ In Fig. ^ we display the 
phase diagram obtained with the two-site approximation 
and compare with exact diagonalization (ED)iiLi^ The 
agreement is remarkably good, in particular concerning 
the value of fj^ 2 - Note that both methods are essen- 
tially based on an approximation for the self-energy by a 
rational function, but the two-site approximation uses a 
function of much simpler structure, i.e., with fewer poles. 

The extension of these results to T > is best de- 
scribed in the grand canonical ensemble. In Fig. [21 four 
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FIG. 2: (U, T, fi) phase diagram for t% — in two-site SFA. 
The plot is constructed from the (U, fi) phase diagram at T = 
(Fig. EJ .the (T, U) phase diagram at half-filling (fi = U/2) 
as in Ref. l3a . and two (T, fj,) sections at fixed interaction U = 
5.3 1* and U = 6£*. Due to particle-hole symmetry the phase 
diagram is symmetric with respect to fi = U /2. Only the 
electron-doped region > U/2) is displayed. 



representative planes in the (U,T,n) phase diagram are 
shown: The (T, U) plane at half-fillingS (fx = U/2), the 
{U,/.i) plane at T — discussed in the last paragraph, 
and two (T,fi) planes at fixed interaction. The latter in- 
teractions are U = 5.3t\ and U = 6£*, which are slightly 
smaller and larger than U C 2, respectively. Our results 
show that the metallic and insulating regions are sepa- 
rated by a first-order transition surface T C (U, fj,) which 
coincides with the transition line T C (U) at /i = U/2 and 
with nf 2 at T = 0. The coexistence region associated 
with this transition lies between two surfaces T cl (£/,//) 
and T C 2(U, /i) where the metastable insulating and metal- 
lic solutions disappear. They intersect the T = plane in 
the lines nf^U) and n%{U) and meet in a line of second- 
order critical endpoints. 

To detect phase separation we calculate the density 
n(n) for given U and various T > as shown in Fig. [3] 
For the first-order filling-controlled MIT at T > we 
find hysteresis, i.e., the densities n™ et and in the two 
phases at the transition \i c arc different. The same behav- 
ior has been obtained in quantum Monte Carlo (QMC) 
DMFT calculations?i2i2i and with the dynamical cluster 
approximation (DCA) for the two-dimensional Hubbard 
modeli^i In the latter work the similarity of the n(/z) 
curves to the p-V isotherms of a van-der-Waals gas was 
noted. The insulating solution is characterized by a low 
charge compressibility and thus corresponds to the (in- 
compressible) liquid, while in the metal x M stays finite at 
the transition even for T — > 0. 

If the system is prepared at a density within the dis- 
continuities [n™ s , n™ 6 *], the free energy is minimized 
by the formation of a phase mixture with a fraction 
%met — (n — n" ls )/(n™ et — n" ls ) of the metal and 1 — x met 



FIG. 3: Density as a function of the chemical potential for 
U = 6tl and ti — 0. The temperature values are given in units 
of 10 -2 t*. Solid lines correspond to the stable phases, while 
dotted lines correspond to the metastable metallic (fi < fj, c ) 
and insulating (fj, > fj, c ) phase, as well as to a third stationary 
point of the grand potential functional, which does not cor- 
respond to any physical phase because Xm < ( m analogy to 
the p-V diagram of the van-der-Waals gas). The discontinuity 
of n at the transition from the larger value in the metal to a 
value close to n — 1 in the insulator is indicated by an arrow. 

of the insulator. With QMC it is difficult to go beyond a 
calculation of n(/i) and determine also fj, c , n™ e4 , and n™ s . 
Within the SFA the latter quantities can easily be calcu- 
lated, since this requires only a comparison of the grand 
potential (per site) f2 in the two phases, a quantity that 
is calculated with high precision anyway. This allows us 
to go from the {T,y) phase diagram to the (T,n) phase 
diagram, which displays the region between n™ e4 (T) and 
n ins (rp-^ w j iere th. e system is unstable against phase sep- 
aration (Fig. Note that phase separation occurs only 
away from half-filling. This will be discussed further in 
the next section. 



IV. NN AND NNN HOPPING 

In the following section we investigate the influence of 
NNN hopping on the paramagnetic phase diagram. We 
focus on the case of t\/t\ = 3/7. Other values of ti yield 
qualitatively similar results. The energy scale is set to the 
square root of the second moment^ of the noninteracting 
DOS, t* = y/tf +tf. 

A. Filling-controlled transition 

For £2 7^ the surface T c (U,fi) of the first-order tran- 
sition in the (U, T, fi) phase diagram retains the overall 
shape as in Fig. [21 but is no longer symmetric with re- 
spect to n = U/2 due to particle-hole asymmetry. The 
value of the critical interaction U C 2 turns out to be rather 
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FIG. 4: (T,n) phase diagram for U = 5.3 1\ and [/ = 6t*; 
both for £2 = 0. Small circles indicate the critical point, 
corresponding to the second-order endpoint of the transition 
line in the (T, fj,) phase diagram. The small dip in n™ s for 
U = 6t* (arrow) corresponds to an intermediate phase and 
might be an artefact of the two-site approximation. 



insensitive to the value of for t^/t* = 3/7 we find 
U c w 5.45i*, in comparison to U c « 5.88^ for pure NN 
hopping. Also the (T, n) phase diagrams look very sim- 
ilar to that for of pure NN hopping (Fig. EJ). However, 
particle-hole asymmetry is immediately evident by the 
very different sizes of the two phase-separated regions, 
which enclose the Mott transition driven by hole- and 
electron-doping, respectively. 

In addition to the phase-separated region, also the 
spinodal curves T c i(n) and T C 2(n) are included in Fig. [S] 
To clarify their meaning let us consider the behavior of 
the system upon heating at given density n = ?io, start- 
ing from the metallic phase at T = (arrow in Fig. |5]i). 
When no = n™ et (T) the system becomes thermodynam- 
ically unstable against the formation of a phase mixture 
containing the metallic phase (n > no) and the insulat- 
ing phase (n w 1). However, the purely metallic phase 
(?i = no) can be superheated up to the spinodal curve 
T c2 . The charge compressibility diverges at the spinodal 
curve, indicating increasing density fluctuations in the 
metastable metallic phase. For T > T c2 a purely metallic 
solution no longer exists. Similarly, the insulating phase 
can be supercooled down to T c \. 

For the smaller value of the two interactions in Fig. 
a closer look at the interval around n = 1 reveals a novel 
situation compared to the case of t 2 = (cf. Fig.|5J). Ob- 
viously, phase separation now occurs at half-filling, but 
disappears at a special density no 7^ 1. To understand 
these numerical results it is convenient to go back to the 
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FIG. 5: (T,n) phase diagram for tj/tj = 3/7. (a) U = 
5.514 4* and (b) U = 5.022 t* . The dotted line is the ther- 
modynamic transition line (n™ e and ric™ s )- Inside this line 
(shaded region) the system is unstable against phase separa- 
tion. The arrow in (a) is referred to in the text. 



(T,/i) representation. For given U and t 2 = 0, the tran- 
sition line T c (a*) has a characteristic U-shape (Fig. |5J), 
where the minimum is at /i = U/2. Precisely at this 
minimum no phase separation is observed. This is clear 
from the Clausius-Clapeyron equation 
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which relates the slope of T c (/j,) to the density n™ et (n* ras ) 
and the entropy (per site) S™ et (Sl ns ) in the two phases 
at the transition. The U-shape has a simple thermody- 
namic explanation: The stable phase is determined by 
its lower grand potential (per site) fi = E — nfj, — TS. 
For U < U c2 we know that the metal is stable at T = 0. 
Let now \x = fxi be such that n met = 1 for T = 0. Wc 
can expand the difference Afl(u) = (Q, met — £T" is )t=o as 



Afi(u) = const. 



(10) 
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FIG. 6: Closeup of Fig. EJd near n = 1, for a discussion see 
text. 



where we used dtl/d/j, = — n as well as n™^ = 1. The 
metal thus becomes more stable as n deviates from n = 1. 
However, in DMFT the insulator has a large entropy even 
at T = since there is no ordering of the magnetic mo- 
ments. This entropy gain stabilizes the insulator at a cer- 
tain temperature T > and causes the first-order transi- 
tion. The minimum of T c directly reflects the minimum 
in Afi(/z). All this remains true for <2 7^ 0. However, 
while the minimum of T c lies at n = 1 as a consequence 
of particle-hole symmetry for t% = 0, it shifts to rip ^ 1 
for t 2 ^ 0. 



B. Half-filling (n = 1) 

Phase separation makes the calculation of the phase 
diagram for a fixed density considerably more difficult, 
ft is impossible to decide if a system at some given den- 
sity prefers to assume an inhomogcneous state by inves- 
tigating only that particular density. Instead, we con- 
struct the {T,U) phase diagram for half-filling from sev- 
eral (T,n) phase diagrams (for various U) as in Fig. [5] 
The first-order line of the particle-hole symmetric case 
is then found to be replaced by a phase-separated re- 
gion as shown in Fig. [7] Within this region, the stable 
phase is a mixture of two slightly-doped phases. Their 
composition cannot be inferred from Fig. alone and 
again requires (T,n) sections. A distinctive feature of 
the phase diagram for pure NN hopping at half-filling is 
the triangular-shaped coexistence region (see Fig- El - By 
contrast, in the presence of frustration we find regions of 
phase coexistence, as well as a small region (close to the 
critical point) in which both metal and insulator phase 
cannot exist as pure phases at n = 1 . This is because the 
spinodal lines cross once and meet at the critical point, 
which is indicated by a circle in Fig0 

For a closer look at the critical region consider Fig. [S] 
When U is lowered the transition shifts to larger tcmpcr- 



FIG. 7: (T, U) phase diagram for half filling and 1%/t* = 3/7. 
For comparison, the corresponding plot for pure NN hopping 
obtained within the two-site approximation is shown in the 
inset. The circle denotes the critical point. Note that the 
spinodal curves cross once within the phase separated region 
(arrow). 

atures and becomes more and more asymmetric (Fig.|SJi), 
until it finally occurs only in the electron-doped regime 
(Fig. |Ht>) . For a given U, the first-order line T c (/x) has 
two critical endpoints. In general the density at these 
points, i.e., the critical density n cr i t (U) is different from 
one. However, when n cr n = 1 (which occurs between 
U = 4.69 1* and U = 4.7 1* in the present case), there 
is a critical point in the (T,U) phase diagram for half- 
filling. This point is then characterized by a diverging 
charge compressibility \n- 

Quantitatively, NNN hopping leads to a considerable 
suppression of the transition region to lower tempera- 
tures, while the critical interaction does not change very 
much. How much this will be modified when one goes 
beyond the two-site approximation remains to be in- 
vestigated. Nevertheless, the fact that NNN hopping 
reduces the Mott transition temperature implies that 
the competition between paramagnetic and antifcrro- 
magnetic phases for t\-t2 hopping is more complicated 
than previously expected. 

V. CONCLUSION 

We employed the self-energy functional approach in 
combination with dynamical mean-field theory to com- 
pute the metal-insulator transition in the (U,T,n) phase 
diagram of correlated electrons. To this end the single- 
band Hubbard model with NN and NNN hopping on a 
Bethe lattice was investigated in the paramagnetic phase. 
For pure NN hopping we find a transition scenario that 
is consistent with previous calculations, namely a first- 
order filling-controlled metal-insulator transition with 
phase separation. No phase separation occurs at half- 
filling, where the transition as a function of the intcrac- 
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(b) U/t*=4.690 



1 1.001 1.002 1.003 1.004 1.005 
n 

FIG. 8: (T,n) phase diagram for 1%/t* = 3/7. (a) U = 4.7i* 
and (b) U = 4.69 1*. 



T] — > + . Close to the critical point the phase diagram 
for half-filling and ti ^ resembles that for £2 = and 
low doping. We note that for the Hubbard model with 
pure NN hopping in d = 2 phase separation was excluded 
analytically at half-filling jsi but this argument does not 
apply to d > 3. 

It is clear that in reality the long-range nature of 
the Coulomb interaction will counteract the formation 
of a strongly inhomogeneous charge distribution. Nev- 
ertheless the instability with respect to phase separation 
found in the Hubbard model hints at intriguing physics. 
Namely, it shows that the inclusion of further interac- 
tions will lead to new and interesting phenomena, e.g., 
more complex ordering phenomena. 

Another open question is the influence of NNN hop- 
ping on the antifcrromagnetic phase. QMC results in- 
dicate the suppression of antiferromagnetism42^ How- 
ever, to decide whether this will be strong enough to 
reveal the metal-insulator transition in the paramagnetic 
phase, the SFA needs to be evaluated for larger reference 
systems. The present results suggest that frustration not 
only suppresses the temperature of the antiferromagnetic 
transition but also that of the Mott transition. Appar- 
ently NNN hopping leads to a more subtle competition 
between metallic, insulating, and magnetically ordered 
phases than previously thought. Therefore it remains a 
challenging task to understand the physical properties of 
strongly correlated materials such as V2O3 in terms of a 
minimal electronic correlation model. 
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